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Analytically tractable dynamical systems exhibiting a whole range of normal and anomalous deterministic 
diffusion are rare. Here we introduce a simple non-chaotic model in terms of an interval exchange transfor¬ 
mation suitably lifted onto the whole real line which preserves distances except at a countable set of points. 
This property, which leads to vanishing Lyapunov exponents, is designed to mimic diffusion in non-chaotic 
polygonal billiards that give rise to normal and anomalous diffusion in a fully deterministic setting. As these 
billiards are typically too complicated to be analyzed from first principles, simplified models are needed to 
identify the minimal ingredients generating the different transport regimes. For our model, which we call the 
sheer map, we calculate all its moments in position analytically under variation of a single control parameter. 
We show that the sheer map exhibits a transition from subdiffusion over normal diffusion to superdiffusion 
under parameter variation. Our results may help to understand the delicate parameter dependence of the 
type of diffusion generated by polygonal billiards. We argue that in different parameter regions the transport 
properties of our simple model match to different classes of known stochastic processes. This may shed light 
on difficulties to match diffusion in polygonal billiards to a single anomalous stochastic process. 


PACS numbers: 05.40.-a, 45.50.-j, 02.50.Ey, 05.45.-a 


Consider equations of motion that generate dispersion of an ensemble of particles as the dynamics evolves 
in time. A fundamental challenge is to develop a theory for predicting the diffusive properties of such a system 
starting from first principles, that is, by analyzing the microscopic deterministic dynamics. Here we introduce 
a seemingly trivial toy model that, analogously to polygonal billiards, exhibits dispersion but is not chaotic 
in terms of exponential sensitivity with respect to initial conditions. We show that our simple map model 
generates a surprisingly non-trivial spectrum of different diffusive properties under parameter variation. 
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Figure 1: Example of a polygonal billiard channel in which single point particles scatter elastically with sawtooth walls [3, 12, 20, 21]. Shown 
is how a beam of particles is split by the corners of the billiard while propagating. 


I. INTRODUCTION 


How macroscopic transport emerges from microscopic equations of motion is a key topic in dynamical system theory and 
nonequilibrium statistical physics [1-7]. While microscopic chaos, characterized by positive Lyapunov exponents, typically 
leads to Brownian motion-like dynamics by reproducing conventional statistical physical transport laws, for weakly chaotic 
dynamical systems where the largest Lyapunov exponent is zero the situation becomes much more complicated [3, 8 , 9]. Such 
non-trivial dynamics is relevant for many topical applications like, for example, nanoporous transport [10-14]. In the former 
case the mean square displacement (MSD) of an ensemble of particles grows linearly in the long time limit, ( x 2 ) ~ f 7 with 
7=1 defining normal diffusion. In the latter case one typically finds anomalous diffusion with 7 ^ 1, where for 7 < 1 one 
speaks of subdiffusion, for 7 > 1 of superdiffusion [ 8 , 9, 14, 15]. 

To our knowledge only a few deterministic dynamical systems are known exhibiting all three regimes of subdiffusion, normal 
diffusion and superdiffusion under parameter variation. Examples of one-dimensional maps are a Pomeau-Manneville like 
model where anomalous diffusion originates from an interplay between different marginally unstable fixed points [16]. The 
climbing sine map displays exactly three different diffusive regimes with 7 = 0 , 1 , 2 corresponding to periodic windows and 
chaotic regions connected to period doubling bifurcations and crises [17, 18]. For the two-dimensional standard map numerical 
evidence exists for a transition from sub- to superdiffusion generated by a mixed phase space [19]. Least understood is diffusion 
in two-dimensional polygonal billiards [3, 8 , 12, 13, 20-22], see Fig. 1 for an example. By definition these systems exhibit 
linear dispersion of nearby trajectories with zero Lyapunov exponents for typical initial conditions, hence non-chaotic behavior. 
However, they nevertheless generate highly non-trivial dynamics due to complicated topologies yielding pseudohyperbolic fixed 
points and pseudointegrability. For this reason they are sometimes called pseudochaotic [3, 8 ]. A line of numerical work on 
periodic polygonal billiard channels revealed sub-, super- and normal diffusion depending on parameter variation [12, 13, 23-26]. 
Rigorous analytical results are so far only available for periodic wind-tree models supporting an extremely delicate dependence 
of diffusive properties on variation of control parameters [27, 28]. For the mathematical derivations it has been exploited 
that polygonal billiards can often be reduced to interval exchange transformations (IETs), see also [29-31], These are one¬ 
dimensional maps generalising circle rotations which cut the original interval into several subintervals by permuting them non- 
chaotically. Both polygonal billiards and IETs are known to exhibit hihgly non-trivial ergodic properties, and in general there 
does not seem to exist any theory to understand the complicated diffusive dynamics of such systems from first principles. 
Random non-overlapping wind-tree models and related maps, on the other hand, enjoy a kind of stochasticity which appears 
analogous to the dynamically generated randomness of chaotic systems, leading to sufficiently rapid decay of correlations and 
good statistical properties. Consequently, these models have been found to yield normal diffusion that is indistinguishable from 
Brownian motion [32, 33]. In contrast, periodic polygonal billiards have very long lasting dynamical correlations and poor 
statistical properties, which are associated with very sensitive dependence of their transport properties on the details of their 
geometry [ 12 ], 

These difficulties to understand diffusion in polygonal billiards on the basis of dynamical systems theory are paralleled by 
difficulties in attempts to approximate their diffusive properties by stochastic theory: There still appears to be a controversy in 
the literature of whether continuous time random walk theory and Levy walks, fractional Fokker-Planck equations or scaling 
arguments should be applied to undestand their anomalous diffusive properties, with different approaches yielding different 
results for the above exponent 7 of the MSD [3, 8 , 21, 34]. While all these theories are based on dynamics generated from 
temporal randomness, spatial randomness leads yet to another important class of stochastic models, called random walks in 
random environments, which yields related types of anomalous diffusion: An important example in one dimension is the Levy 
Lorentz gas where the scatterers are randomly distributed according to a Levy-stable probability distribution of the scatterer 
positions. This model has been studied both numerically and analytically revealing a highly non-trivial superdiffusive dynamics 
that depends in an intricate way on initial conditions and the type of averaging [35-37]. This work is related to experiments on 
Levy glasses where similar behavior has been observed [38]. 
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Figure 2: Space-time plot illustrating the action of the one-dimensional slicer map S a for a = 1 /3 defined by Eqs. (1,2), where m (horizontal 
axis) denotes space and n (vertical axis) time: Shown is the diffusive spreading of points that at n = 0 are uniformly distributed on the unit 
interval centered around m = 0. As the map is one-dimensional the columns are only a guide to the eye. This map is designed to mimic the 
mechanism of beam-splitting in polygonal billiards depicted in Fig. 1. 


Motivated by the problem of understanding diffusion in polygonal billiards, in this paper we propose a seemingly trivial non- 
chaotic map by which we attempt to mimic the dynamics illustrated in Fig. 1 : Shown is a beam of point particles and how it splits 
due to the collisions at the singularities (corners) of a polygonal billiard channel. This mechanism is intimately related to the 
connection between polygonal billiards and IETs referred to above. We thus try to capture this slicing dynamics by introducing 
a specific IET defined on a one-dimensional lattice, see Fig. 2. Here the loss of particles propagating further in one direction is 
modeled by introducing a deterministic rule following a power law for the jumps from unit cell to unit cell. 

This simple non-chaotic model, which we call the slicer map, generates a surprisingly rich spectrum of diffusive dynamics 
under parameter variation that includes all the different diffusion types mentioned above. We mention in passing that it provides 
another example where normal diffusion is obtained from non-chaotic dynamics. However, differently from the cases of [32, 33] 
and analogously to periodic polygonal billiards, it is completely free of randomness. Our simple model might help to understand 
why the type of diffusion in polygonal billiards is so sensitive under parameter variation. It might also shed some light on the 
origin of the difficulty to model polygonal billiard dynamics as a simple stochastic process. 

Our paper is organized as follows: In Section 2 we define the slicer model and analytically calculate its diffusive properties. 
To the end of this section we study our model analytically and illustrate it numerically for a parameter value that is characteristic 
for the dynamics. In Section 3 we compare the deterministic slicer dynamics with existing stochastic models of anomalous 
diffusion. Section 4 contains concluding remarks. 


II. THE SLICER DYNAMICS 
A. Theory 

Consider the unit interval M := [0,1], the chain of such intervals M := M x Z, and the product measure fi := A x 8z on M, 
where A is the Lebesgue measure on M and 8z is the Dirac measure on the integers. Denote by ttm and 7Tz the projections of 
M on its first and second factors. Let x be a point in M, X = ( x , m ) a point in M, and M m := M x {to} the m-th cell of M. 
Subdivide each M m in four sub-intervals, separated by three points called “slicers”, 

{1/2} x {to} , {£ m } x {to} , {1 - l m j x {to} , 

where 0 < i m < 1/2 for every to £ Z. 

The slicer model is the dynamical system (M, p, S) which, at each time step n £ N, moves all sub-intervals from their cells 
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to neighbouring cells, implementing the rule S : M —>• M defined by 


S(x, to) 


(cc, to — 1) if 0 < x < £ m or i < x < 1 — £ m , 
(x, to + 1) if < x < ^ or 1 — i m < a; < 1. 


For every a > 0, let us introduce the family of slicers 


L a 



1 -^j(a)) 




1 

(|j| + 2 1 /“) a ’ 



(1) 


( 2 ) 


The sheer map is denoted by S Q if all slicers of Eq.(l) belong to L a : t rn = £ m (a). Obviously, for every a > 0, S a preserves 
fi and is not chaotic: Its Lyaponuv exponent vanishes, as different points in M neither converge nor diverge from each other in 
time, except when separated by a sheer in which case their distance jumps. This is like for two particles in a polygonal billiard 
where one of them hits a corner of the polygon while the other continues its free flight, see Fig. 1. But the separation points 
constitute a set of zero ft measure, hence they do not produce positive Lyapunov exponents. 

The dependence of the dynamical rule Eq.(l) on the coarse grained position in space to is a crucial aspect of the sheer model, 
which distinguishes it from ordinary IETs. That the slicers get closer and closer to the boundaries of the cells when the absolute 
value of m grows is meant to reproduce, in a one-dimensional setting, what is illustrated in Fig. 1: The corners of periodic 
polygonal billiards split beams of particles into thinner and thinner beams as they travel further and further away from their 
initial cell. In two dimensions this operation is fostered by the rotations of the beams of particles, something that is not possible 
in a single dimension. This thinning mechanism due to slicing is mimicked by the power law dependence in Eq.(2). In effect 
this means that our sheer particles perform a deterministic walk in a Levy potential. This quite trivial setting has, as we shall 
see, rather non-trivial consequences. The power law dependence is a mere assumption at this point in order to define our model. 
It would have to be developed further to move the sheer map closer to actual polygonal billiard dynamics. 

The diffusive properties of the sheer dynamics will be examined by taking an ensemble of points Eq in the central cell 
Mo = M x {0} and studying the way S a spreads them in M. One finds that in n time steps the points of Eq reach M n and 
M_„, and that the cells occupied at time n have odd index if n is odd, and have even index if n is even. 

More precisely, taking 


P n = {j £ Z : j is even and \j\ < n}, D n = {j £ Z : j is odd and \j\ < n} 


we have 

5” M 0 = [J [Rj x {j}) if n is even , S™ M 0 = [J (i?. 7 x {j}), if n is odd 

j&Pn j£D n 

where Rj x {j\ C Mj, and Rj C M is an interval or a union of intervals if Eq = Mq, with Ri 0 Rj = 0 if i ^ j. 
Let dv o := po(X)dfi be a probability measure on M with density 


Po 



1, 

0, 


ifX £ Mo 

otherwise . 


(3) 

(4) 


(5) 


This measure evolves under the action of S a describing the transport of an ensemble of points initially uniformly distributed in 
Mo- In the following we will always adopt this initial setting, which mimics a ^-function like initial condition as is common in 
standard diffusion theory, adapted to a lattice by filling a unit cell in it. If the initial condition were confined within M m with 
to ^ 0, nothing would change qualitatively. However, if we would fill a unit cell non-uniformly with particles, e.g., by choosing 
points close to the boundary of a cell, clearly we would observe very different dynamics. Hence, there is dependence of the 
outcome on the initial measure as is typical for IETs. Here we characterize its dynamics by choosing a sufficiently ‘nice’ initial 
measure. 

Requiring conservation of probability, the evolution //,, of uq 
R C M. Its density is given by 

Pn{X) = | J 

In the line above Eq.(6) S~ n is intended in the set-theoretical 
restricting to the specific initial condition given by cell M 0 , and to the part of M that the points initially in M 0 reach at any 
finite time n, the preimage of a point is a single point and the inverse of the map is defined as follows: Consider the evolution 


at time n is given by v n (R) = vo(S a n R) for every measurable 

ifXGS^M 0 (6) 

otherwise . 

sense, since S~ l X is not a single point, in general. However, 
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of the initial distribution N^ = S™M 0l n = 0,1, 2,3,..and define the maps T n = S a |j^( n) : —> N ( ra+1 ) (we drop 

a from T n for the sake of notation). These maps are surjective. They are also injective. Indeed, suppose x\ , x% G N( n ) yields 
T n {x i) = T n (x 2 ), then ttm(T„(xi)) = ir m(T u (x 2 )) =: £ and, since S a does not change the first component of any point x 
(i.e. 7 tm{S(x)) = 7 we have that also ttm{%i) = tvm{x 2 ) = £■ But in there is only one point with first component 
£ for any £ £ (0,1), thus x\ = X 2 ■ The map that gives the evolved distribution at time n is 

© a,n ■= T n -1 o T„_2 o • • • o Tq : Mq —» M, ( 7 ) 


and the dynamics is given by the family of invertible maps {© a ,n}neN- Obviously © a ^ = T 0 1 o • • • o T n _ 2 o T n \. Since 
TYM{T n (x)) = 7 tm{x), for any A C N n we have A(7 tm{T u (A))) = A(7Tm(^ 4)), where A is the Lebesgue measure. In the same 
way if A C N^ n+1 \ then X(ttm(T~ 1 (A))) = X(ttm(A)). From this it follows that if A C N^°\ then A(7TM(©n(^))) = 
X(ttm(A)) and if A C then X(ttm (©aii(^))) = X(ttm(A)). In other words, maps {© a ,n}raGN preserve the Lebesgue 
measure and jl is also invariant w.r.t. the same family of maps. 

Apart from being formally precise with defining the inverse for our model, it is interesting to conclude that S' -1 depends on 
the initial condition. More physically speaking, this appears to be a consequence of the spatial translational symmetry breaking 
with respect to the different slicer positions in the different cells. The lack of a more general definition of S’ -1 implies in turn 
that the slicer map is not time reversible invariant in the sense of the existence of an involution [1-3], however, we don’t need 
the latter property for our calculations. 

Consider now the sets 

Rj := S'" M 0 n Mj = © Q ,n(M 0 ) n M), j = -n, ( 8 ) 

which constitute the total phase space volume occupied at time n in cell Mj. Their measure 

Aj := fi(Rj) = X(TT M (Rj))Sz(j) = X(n M {Rj)) (9) 

equals the probability u n (M :j ) of cell j at time n : as /i is {© ajrl }„ e j!j-invariant and ©«„ are invertible, we have 

Aj = KRj) = iXG-^Rj)) = A(Mo n e^ n (M s )) = ^(©"U^)) = = Mm , 

and E"=-n A?' = A( u ”=-n©a,n( M o) H Mj) = = jx{M 0 ) = 1. Indeed, 6“U M o) n Mj = 0 for \j\ > n 

and U^L_ 00 Mj = M. In other words, the A,’s define a probability distribution which coincides with v n {Kj, 1 ) and, thus, is a 

marginal probability distribution of v n . Starting from the “microscopic” distribution v n on M, we can now introduce its coarse 
grained version f) ( ;j as the following measure on the integer numbers Z: For every time n G N, the coarse grained distribution is 
defined by 


n G (A = I Aj if 7 e 

rn \J / | q otherwise . 


( 10 ) 


A_ n and A n are called traveling areas , Aj is called sub-traveling area if \j\ < n. 

Remark 1. From the definition of S a and the initial condition Eq.(5), we have Aj = A—j for all j G Z. Thus, p G (j) is even, 
PnU) = Pn ( — 7 )> an d i fs odd moments vanish. 

The coarse grained distribution will be used to describe the transport properties of the coarse grained trajectories 
{7Tz(S'”A’o)} n gN C Z, with Xq G Mq. This way, p G becomes the discrete analog of the mass concentration used in ordi¬ 
nary and generalized diffusion equations for systems that are continuous in time and space [9, 15]. Accordingly, we can define 
a discrete version of the MSD as 


(A X 2 n ) := J2 A 0 j 2 (11) 

j=—n 

for p G , where j is the distance travelled by a point in Mj at time n. Then, for 7 G [0, 2] let 

Ml) := lim . (12) 

n—> 00 Tl* 

If T a { 7*) G (0, 00 ) for 7* G [0, 2], 7* is called the transport exponent of the slicer dynamics, and T a ( 7 *) yields the generalized 
diffusion coefficient [9, 15]. 
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Remark 2. Due to symmetry the mean displacement (AX n ) := Y^j=-n AjJ vanishes at all n, hence there is no drift in the 
slicer dynamics. 

Note that A 3 equals the width of the interval Rj, which is determined by the position of the slicers in the y-th cell, once a is 
given. Therefore, A :j can be computed directly from Eq.(2). For the traveling areas we have 

An=£n ~ 1 = (\n\-l + 2i/«) =A ~ n (13) 

while for the non vanishing sub-traveling areas we have 

Aj = £,3 ‘ hl ~ iul+1 = (|j| — 1 + 2 1 / Q ) a “ (|j| + l + 2Va)“ ’ (14) 

For even n > 2 this implies 


II 

eTs 

Cl 

( 2(4-4) 

) i\2k-l\ ~ ^| 2 fc+l| 

1 ^\n—l\ 

\ 0 

, for j = 0 

, for \j\ = 2k, k= 1 ,.. 

, for jjj = n 
, elsewhere 

n -2 
• ’ 2 

while for odd n > 3 it implies 




II 

Oe 

CL 

^| 2 fe| — ^| 2 fc+ 2 | , 
i-|n— 1 | ■ 

for bl = 2fc T 1, k = 0,.. 
for \j\=n 
elsewhere . 

n —3 
. . , 2 


(15) 


(16) 


Remark 3. Using Eq.(2) in Eqs.(15,16) for large n, one obtains that the tail of the distribution (large j) goes (independently of 
the parity of n) like jfrf (j) ~ 2a/|j| a+ 1 I/u-i <n }, i.e. p ® has heavy tails. Note that for a G [0, 2) these tails correspond exactly 
to the ones of a Levy stable distribution [9], However, for j = ±n the probability is much larger, p'f (n) ~ l/|n|“. 

We are now prepared to prove the following result: 

Proposition 1. Given a G (0,2) and the uniform initial distribution in Mq, we have 


T a ( 7) 


Too if 0 < 7 < 2 — a 
5 ^ if y = 2 — a 
0 if 2 — a < j < 2 


(17) 


hence the transport exponent takes the value 2 — a with (AX 2 ) ~ n 2 a . For a = 2 the transport regime is logarithmically 
diffusive, i.e. 

(AX 2 )~\ogn (18) 


asymptotically in n. 

Proof. Because of the symmetry of p„ let us consider only the cells Mj with j G N, so to obtain: 


n -j / n —1 

T a {l) = 2 lim — AD 2 = 2 lim — ( V'' AH 2 
V " n-»oo n't ^ 3 moont \ ^ 3 

3=0 \j =0 


■ A r 


Because of Eq.(13), the travelling area yields 


lim 


n—too (n 4- 2 1 / a — l) c 


oo if 0 < 7 < 2 — a 
— = ^ 1 if 7 = 2 — a 

nl 1 0 if2 — a<7<2 


For the sub-travelling areas, by introducing Q n := Ajj 2 , we will show below that 


i • Qn 

n—H>o n't 


oo if 0 < 7 < 2 — a 


a 


if 7 = 2 — a 


2 — a 

0 if 2 — a< 7<2 


(19) 


( 20 ) 


( 21 ) 
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Remark 4 . Note that the traveling and the sub-traveling areas produce exactly the same scaling for the MSD. We will come 
back to this fact in Section III. This result can also be obtained by calculating the second moment of the probability distributions 
of these two different areas directly from the expressions given in Remark 3 above. 

To prove Eq.(21), consider that the series Q n assumes a different form depending on wheter n is even or odd. If it is even and 
larger than 2 we have 


n— l ir -1 

Qn = Ajj = 4 

j=o,jeP n j=i 


( 2 j + 2 1 /" — l) ( 2 j + 2 V“ + l) 


This sum has a telescopic structure that allows us to rewrite it as 


£-i 


Qn = 4 X] 


2j-l 


(n — 2) 2 


^ (2 j - 1 + 2 V«)“ (n - 1 + 2 V“) c 


( 22 ) 


(23) 


Let R n be the first term of Q n . Introducing f(j) := 


2.7-1 


(2j-l+ 2 i/«)“ 


we can write 


f? 


n 



2j -1 

(2j - 1 + 2!/“)“ 


4 E m 

3 =1 


(24) 


The derivative 


= 2 [2(1 — a)j + 2 1 / a + a — l] 
■' (2j + 2V“ - 1)“ +1 


(25) 


shows that / is increasing for 0 < a < 1, while for 1 < a < 2, / grows for j < j(a) and decreases for j > j(a), with 
j(a) = (1 — a — 2 1//a )/2(l — a). For 0 < a < 1 / is strictly increasing, hence 


r* —— 1 
r o J- 


f -1 


f{x)dx < E /O') < / /OOdx . 

■ /u i=i Jl 

We have to distinguish two cases, a < 1 and ct = 1. In the first case we have 


f(x)dx = | 


(n — 3 + 2 1 / a ) 2 ~ a _ 1/a _ (n-3 + 2 1 / a ) 1 - 


(V /a _ ^ 


2 - a 
i-“ 2 1 /" - a 4-1 


1 — a 


and 


f(x)dx = i 


(2 — a)(l — a) _ 

(n - 1 + 2 1 /“) 2_ “ _ 1/a _(n - 1 + 2 1 /“)i- 


+ (2 i/a + 1) 

therefore taking the n —> oo limit we have 


2 — a 

i/a , ,m-«. (2 17 “ + a-l) 


1 — a 


+ 


(2 — a)(l — a) 


1 , t -± | oo if 0 < 7<2 — a 

lim — / f(x)dx = { 2 ( 2^1 if 7 = 2 - a 
n^oon'TJo > 


and 


1 f 7 

lim — f(x)dx = 

n—toc n~ Ji 


0 if2 — a<7<2 


oo if 0 < 7 < 2 — a 
2(2^y if 7 = 2 — a 

0 if 2 -ct< 7<2 . 


(26) 


(27) 


(28) 


(29) 


(30) 


7 























For a = 1 the two integrals differ, but the bounding limits coincide. Therefore one obtains 


,. Rn 
hm — = 

n—> oo 77,7 


lim — 

n—>oo 77,7 


E m 

j=i 


oo if 0 < 7 < 2 — a 

- if 7 = 2 — a 

2 — a 

0 if 2 — a< 7<2 


( 31 ) 


If 1 < a < 2, f decreases for j > j(a). Hence, introducing j a = \_j(a)\ , where is the integer part of x, R n can be 
expressed as: 


Rn =4 


9-i 




o=i 


j=j C+l 


(32) 


Dividing by n 7 and taking the n —> oo limit, the first term vanishes for all 7 > 0 while the second term can be treated as above 
to obtain the same as Eq.(31). Recalling Eq.(23), this eventually implies Eq.(21). For odd n, one proceeds similarly. 

In summary, the MSD grows like (AX 2 ) ~ n 2 n for a £ (0,2), and the trivial sheer map S a enjoys all power law regimes 
of normal and anomalous diffusion as a varies in ( 0 , 2 ). 

Remark 5. From Eq.(2) it follows trivially that for ct-tOwe have tj = 1/2, Vj £ Z. This means that everywhere on the slicer 
lattice half unit intervals are mapped onto half unit intervals in neighbouring cells in the same direction of the previous jump 
generating purely ballistic motion. Consequently for a = 0 the MSD grows like (AXf) ~ n 2 {n —> 00 ). 

Repeating the previous reasonings by computing the correspondingly different integrals for a = 2 for a uniform initial 
distribution in Mq we find that 


2M7) 


+00 if 7 = 0 
0 if 0 < 7 < 2 


(33) 


The upper and lower bounds of the integrals corresponding to Eq.(26) feature leading logarithmic terms, which yields Eq.(18). 

□ 


Remark 6. By an analogous calculation, or alternatively by looking at the second moment of the probability distributions, cf. 
Remark 4 above, one can see that for a > 2 

{AX%) const, (n — > 00 ) . (34) 

That is, in terms of the MSD localisation sets in, although from the definition of the slicer map it is intuitively not clear why this 
should happen. 

These results can now be generalised by calculating the asymptotic behavior of the higher moments AX? of 

n 

{AX?) = ]T Ajj* . (35) 

j=-n 

Theorem 2. For a £ (0, 2] the moments (AXjf) with p > 2 even and initial condition uniform in Mq have the asymptotic 
behavior 

(AX?) ~ nP~ a (36) 


while the odd moments (p = 1, 3,...) vanish. 
Proof. We want to compute the limit 


1 1 _ ^ 

L{a,p) := lim —(AX%) = 1™ — E A ^ P ' (37) 

n—too Ji 1 n—>oo 77, / z ' 

j=-n 

As observed in Remark 2, the symmetry of p ( f implies that the sums with odd p to vanish. For the even moments it suffices to 
consider the positive j’s, 


L(a,p) 


lim — 

n—> oo 77,7 


(^T Ajj p + A n rA 


(38) 
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We now prove that 


n— 1 


i A 

lim — V Aif = 

n—t oo 77,7 ^ 


a 


3=0 


p — a 


oo if 0 < 7 < p — a 
>0 if 7 = p — a 
0 if 7 > p — a . 


In order to do so, for even n let us introduce 

n— 1 1 

£ Ajf = 2^ 


j=o,jeP n 

A simple induction procedure leads to 

4-2 


1=1 


(2j + 2 1 /° t — l) (2j + 2 1 /« + l) c 


V = 2 P ■ V + 1 ' )P ~ jP _ 

U (2l + l + 2 1 /«)“ („-l + 2Va) 


(»-2) p _ R («-2f 

a 


(n - 1 + 2 !/a) c 


which defines //,, in terms of addends of the form 

(J + 1) P -J P 


/(j) = 


= £ 


p ta j p - fe 


(2j + l + 2 1 /«)“ Yi\kJ (2j+ l+2 1/a ) 


with derivatives given by 


fe=l V 7 


>\ [2 (p -k- a)j + {p- k)( 1 + 2 1 /“)] fe _ 1 


(2j + 1 + 2 1 /«)«+i 


= £/fc(j) , 


fe=l 


(39) 


(40) 


(41) 


(42) 


(43) 


where 

f n = M [2(p - fc - ct)j + (p - fc)(l + 2 1 / Q )] , p _ fe _ 1 
JfcU ) ( 2j + 1 + 2 1 /“)“+ 1 J 


(44) 


For 0 < a < 1 and all j > 0 we have fk{j) > 0 for k = 1,... ,p — 1, while f p (j) < 0. Because \f p (j)\ < f 

is positive and / increases for all j > 0. For 1 < a < 2 and p = 3 , one has /(j) = (3 j 2 + 3 j + 1)/(2 j + 1 + 2 1 / a )°‘, 
which is increasing for j > 0, while for p > 4 one obtains fk(j) > 0 for k = 1,... ,p — 2, and f p -i(j), f p (j) < 0. Because 
1/p-iGO + f P (j)\ < fi(j) +/2O'). fU) is increasing for j > 0, even for 1 < a < 2. 

Therefore our sum is bounded from above and below by 


1-2 


1-2 


I f{x)dx < Y f(j) < / 
Jo j=0 J 1 

for all a £ (0, 2]. Taking the limit as done previously, we eventually obtain 

00 if 0 < 7 < p — a 
if 7 = p — a 


l-i 


f{x)da 


lim = 
n—>■ 00 77 ,7 


P 


2P{p — a) 


lim £ 

n—> oo 77,7 


0 if 7 > p — a , 


for all a £ (0,2]. If n is odd one proceeds similarly to obtain the same result. 
For the travelling area Eq.(13) we have 


00 if 0 < 7 < p — a 

v 

if 7 = p — a 


p — a 

0 if 7 > p — a 


lim A n n p = lim -—-—- 

n —¥oo n—>00 (ji 2 1 / a — l) c 


^ iOoif0<7<p — a 
— = < 1 if 7 =p-a 

n 1 0 if 7 > p — a . 


(45) 


(46) 


(47) 


Hence, the same scaling for traveling and sub-traveling regions as pointed out for the second moment in Remark 4 holds for all 
higher moments. We thus conclude that 


! oo if 0 < 7 < p — a 

p 

- if 7 = p-a 

p — a 

0 if 7 > p — a 

so that the large n behavior of the even moments is given by (AX P ) ~ n p ~ a . 


(48) 

□ 
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B. Example 


In this subsection we illustrate the diffusive transport generated by the slicer map S a for a representative value of a. For this 
purpose we plot the probability distribution using our exact analytical results and compare it to an asymptotic approximation. 
We then draw cross-links to diffusive transport in a polygonal channel. 

As an example, let us consider the case a = 1/3. Why we choose this particular values is explained further below. Here 
we have £,-(1/3) = l/(|j| + 8) 1 / 3 , and the asymptotic behavior of the MSD is given by (AX%) ~ n 5 / 3 , cf. Proposition 1. 
This means that Si / :i is superdiffusive with 7 * = 5/3 and generalized diffusion coefficient X 1/3 = 12/5. From Theorem 2 the 
moments of S // 3 higher than the second have the behavior 

(A X») ~ n p ~ 1/3 . (49) 


The coarse grained distribution for even n reads, see Eqs.(15,16), 


Pn M 


1 

V2k+7 


1 

2 


1 



v 7 3E+9 ! 



0 , 


for m = 0 

for to = 2 k, k = 2 ,..., 
for m = n 
otherwise , 


while for odd n we have 



1 

y2k+8 


v^fc+lO ’ 


for m = 2k + 1 , k = 2 ,..., , 

for m = n 
otherwise . 


(50) 


(51) 


Figure 3 shows the marginal probability distribution function (to) at fixed even n for to > 0, including the last value p 1 ^ (n) 
which is much larger than the values for to close to n. The negative branch of the distribution can be recovered by symmetry. 
Because asymptotically goes like 


P“M 


c 

(to + 2 1 / q )“ +1 ’ 

0 , 


m < n, 
to > n, 


(52) 


where C a is a normalization constant. Fig. 3 compares the numerical values of p c -[ with our asymptotic result for pl/ 3 (m) and 
C 1/3 = 1. Apart from the peak at p ( ;[(n) due to the traveling area, which is covered by Eqs. (50,51), the asymptotic behavior of 
both results is clearly the same. Note that the spike at to = n is analogous to the one found in [36, 37], 

The choice of a = 1/3 is motivated by results on diffusion in the sawtooth polygonal channel studied in [12]. The channel 
geometry is analogous to the one shown in Fig. 1 except that there are no flat wall sections between any two triangles. The 
angle between one side of a triangle and the wall base line has been chosen to 7 r/ 4 . Simulations for this particular case yielded 
a transport exponent of about 7 1 = 5/3, cf. Table 1 in [12], This result is surprising, as naively one would have expected 
irrational polygons to generate transport close to diffusion, and rational polygons to exhibit transport close to ballistic. From 
this perspective the case of 7 r /4 angles should have been perfectly ballistic, while it turned out to be substantially slower than all 
other irrrational cases with parallel walls reported in [12]. This suggests that the main mechanism generating diffusion in this 
channel may have less to do with whether respective polygons are rational or irrational but rather how precisely they slice a beam 
of diffusing particles as modelled by our slicer dynamics. Fig. 2 indicates that in this case diffusion may be slowed down due 
to an incresing fraction of orbits being localized by not contributing to diffusion. This reminds of similar findings for polygonal 
billiard channels presented in [26]. 

We remark that results completely analogous to Fig. 3 are obtained for any other value of 0 < a < 2. This implies that for 
a = 1 our system generates a very strange type of normal diffusion with a non-Gaussian probability distribution. For 1 < a < 2 
it is furthermore surprising that ballistic peaks representing traveling regions are present while the model as a whole exhibits 
subdiffusion. We are not aware of results in the literature where subdiffusive dynamics with coexisting traveling regions has 
been observed. 


III. A SIMPLE STOCHASTIC MODEL OF SLICER DIFFUSION? 

Since deterministic dynamical systems often generate a type of randomness, it is frequently attempted to match their dynamics 
to simple stochastic processes for understanding their transport properties [1-5]. However, as we pointed out in the introduction, 
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Figure 3: Log-log pot of the marginal probability distribution for the sheer map S 1/3 as a function of the position of the m-th cell at fixed 
time n = 10 s . Pn{m) denotes the coarse grained distribution obtained from our exact analytical results Eqs.(50,51) (continuous line). It is 
compared with the asymptotic analytical approximation p“(m), Eq.(52) (dotted line). Apart from the values at m = n, where p„ has a spike 
due to the travelling area, the asymptotic functional forms coincide. 


for diffusion in polygonal billiards such a stochastic modeling turned out to be surprisingly non-trivial [3, 8, 21, 22]. Motivated 
by this line of research, in this section we relate the sheer diffusion to known stochastic models of anomalous diffusion. 

We first summarize our main results for diffusive transport generated by our non-chaotic sheer map under variation of the 
control parameter 0 < a < 2. In the limit of n —> 00 we have: 

1. a = 0: ballistic motion with MSD {AX 2 ) ~ n 2 

2. 0 < a < 1: superdiffusion with MSD {AX 2 ) ~ n 2- “ 

3. a = 1: normal diffusion with linear MSD {AX 2 ) ~ n 

4. 1 < a < 2: subdiffusion with MSD {AX 2 ) ~ n 2 ~ a 

5. a = 2: logarithmic subdiffusion with MSD (AX 2 ) ~ logn 

6. a > 2: localisation in the MSD with {AX 2 ) ~ const. 

Additionally, Theorem 2 gives information about the asymptotic behavior of all higher order even moments scaling as 
{AX P ) ~ n p ~ a in the long time limit forp > 2 and 0 < a < 2. 

As recently highlighted in Refs. [39, 40], there do not seem to exist too many stochastic models exhibiting a transition from 
subdiffusion over normal diffusion to superdiffusion under parameter variation: We are aware of a specific continuous time 
random walk (CTRW) model [16], (generalized) elephant walks [39, 40], and generalized Langevin equations (gLe) including 
fractional Brownian motion [41 — 43]. For these models one can easily check that there is no simple matching between their 
diffusive properties and the above scenario representing sheer diffusion. That is, the scaling of the MSD with parameters by 
switching between all diffusive regimes is generically different from the slicer dynamics, and/or more than one control parameter 
is needed to change the diffusive properties. However, for both the CTRW and the gLe there is a partial matching to the slicer 
diffusion in specific diffusive regimes to which we come back below. Meaning so far we are not aware of any stochastic model 
that (asymptotically) reproduces the slicer diffusion by exhibiting all the different diffusive regimes listed above under single 
parameter variation. 

For all the stochastic models just mentioned the dynamics is generated by temporal randomness, that is, random variables are 
drawn in time from given probability distributions. A second fundamental class of stochastic models is defined by spatial ran¬ 
domness of the positions of scatterers with a point particle moving between them. An important example is the one-dimensional 
stochastic Levy Lorentz gas (LLg) [35-37]: Here a point particle moves ballistically between static point scatterers arranged on 
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a line from which it is either transmitted or reflected with probability 1/2. The distance r between two consecutive scatterers is 
a random variable drawn independently and identically distributed from a Levy distribution, 


A(r) = /3r 


0 1 

0 r /3+l ’ 


r e [ro, +oo) , 


(53) 


where /3 > 0 and ro is a cutoff fixing the characteristic length scale of the system. The LLg shares a basic similarity with the 
sheer map in that its scatterers are positioned according to the same asymptotic functional form as the slicers, see Eq.(2). On 
the other hand, the slicer positions are deterministic while the LLg scatterers are distributed randomly. In more detail, the slicers 
amount to transition probabilities for nearest neighbour jumps on a periodic lattice that decrease by a power law while in the LLg 
the jumps follow a power law distribution with trivial transition probabilites at random scatterers. Finally, the slicer dynamics is 
discrete in time while the LLg is a time-continuous system. 

From these facts it follows that there exists an intricate dependence on initial conditions in the LLg that is not present that way 
in the slicer map: The LLg diffusive properties depend on whether a walker starts anywhere on the line, which means typically 
between two scatterers, called equilibrium initial condition, or exactly at a scatterer, called nonequilibrium initial condition 
[36, 37]. In [36] bounds for the MSD have been calculated in both cases. Interestingly, the lower bound for equilibrium initial 
conditions was obtained from CTRW theory on which we will elaborate further below. 

The results for nonequilibrium conditions have been improved in [37] based on simplifying assumptions by which asymptotic 
results for the position probability distribution of the moving particle could be calculated. It was shown that the LLg only 
displays superdiffusion, which as in the slicer map is governed by traveling (in [37] called subleading) and sub-traveling (in [37] 
called leading) contributions. However, while in the slicer map both these contributions scale in the same way for the MSD, cf. 
Remark 4, for the LLg they yield different scaling laws depending on the order of the moment p and the control parameter /?. 
These different regimes are deeply rooted in the different physics of the system featuring an intricate interplay between different 
length scales. As a consequence, the LLg MSD displays three different regimes with an exponent determined by three different 
functions of /3, cf. Eq.(3) in [37]. All higher order moments could also be calculated for the LLg. Very interestingly, Eq.(13) in 
[37] and all even moments of the slicer map, cf. Theorem 2, exactly coincide by a piecewise transformation between /? £ (0,3/2] 
and a £ (0,1]. In other words, for every /3 in [37] it suffices to fix the slicer’s parameter a so that one of the moments (e.g. the 
second) matches, to match all other moments as well [46]. Then the asymptotic form of the moments ( r p (t )) for all p > 0 is 
given by 


<■ r p it)) 


, if/3 < 1, p < /3 

p(i+p)~p 2 

t 1+3 , if /3 < 1, p > [3 

, if /3 > 1, p < 2/3 — 1 
ti+P-P , if/3 > 1, p> 2/3— 1 


Surprisingly they can be matched with the slicer moments in Eq.(36) by taking 


(54) 


a = 


/3 2 /(l+/3) ifO</?<l 
/3 — 1/2 if 1 < /3 < | 
1 if/3 > | 


(55) 


This means that by using the above transformation both processes are asymptotically indistinguishable from the viewpoint 
of these moments, meaning the slicer map generates a kind of LLg-type walk in the superdiffusive regime if the available 
information on the system (the observables) include the moments only. On the other hand, the transformation is piecewise 
which reflects the different functional forms for the exponent of the moments in the LLg while for the slicer map only one such 
functional form exists. Indeed, it is well known that the moments carry only partial information on the properties of (anomalous) 
transport phenomena, and that knowledge of correlations is necessary to distinguish one class of stochastic processes from 
another [9, 14]. 

Another interesting fact within this context is that for the traveling region alone (called ballistic contribution in [36]) the MSD 
of the LLg scales as ~ t 2- ^ in continuous time t, as was shown in [36]. Formally, this result matches exactly to the MSD of the 
slicer map of ~ n 2 ~ a as calculated in Proposition 1. Note also that the slicer positions generate a probability distribution for the 
sub-traveling region of p^ (j) ~ |/| - “ -1 , see Remark 2, which matches to A(r) of the LLg Eq. (53). 

We now comment on similarities and differences of the slicer diffusion with CTRW theory. For equilibrium initial conditions 
in the LLg it was shown that for 1 < /3 < 2 the MSD is bounded from below by ~ t 3 ~^ [36]. However, this is exactly the result 
for a Levy walk modeled by CTRW theory [16]. Even more, results for all higher Levy walk CTRW moments were recently 
calculated to ~ t p+1 ~^ for p > /3, see Eq.(18) in [44]. With /3 = 1 + a the slicer superdiffusion thus formally (also) matches 
to Levy walk diffusion defined by CTRW theory. On the other hand, from a conceptual point of view a CTRW is constructed 
very differently from both the LLg and the slicer map dynamics. Hence, it is not very clear why a CTRW mechanism should 
apply in both these cases. Another remark is that for the superdiffusive regime of the slicer diffusion of 0 < a < 1 the frozen 
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Levy distribution according to which the transition probabilities have been defined does not belong to the parameter regime for 
which such distributions are stochastically stable in the sense of a generalized central limit theorem [9, 15], which holds only for 
1 < a < 3. This points again at a crucial difference between sheer dynamics and CTRW theory, where for the latter the resulting 
probability distributions are stochastically stable. Our discussion suggests that there is a more intricate interplay between the 
Levy potential we want to mimic, the dynamics we use to obtain it, and the power law distributions we generate from it. 

We conclude this section by a remark on a curious similarity between the sheer diffusion in the subdiffusive regime and 
Gaussian stochastic processes. It was shown in Refs. [41, 42, 45] that for a gLe with power law memory kernels for friction 
and/or noise the MSD exhibits a transition from a constant over ~ log n to subdiffusion in the long time limit. Especially, for an 
overdamped gLe with Gaussian noise governed by a power law anti-persistent memory kernel ~ — t “ 7 the MSD was calculated 
to ~ f 2-7 for 1 < 7 < 2, ~ log t for 7 = 2 and ~ const, for 7 > 2 [45]. Formally these results match exactly to the sheer 
MSD for 1 < a, cf. Proposition 1. However, this overdamped gLe does not exihibit any superdiffusion. And as the probability 
distributions generated by such a gLe are strictly Gaussian in the long time limit there is, again, a clear conceptual mismatch 
to the sheer dynamics: In the gLe the subdiffusion is generated from power law memory in the random noise (by calculating 
the MSD via the Taylor-Green-Kubo formula [41, 42, 45]) while for the sheer map the anomalous dynamics was calculated 
from non-Gaussian probability distributions. We remark that so far nothing is known about the correlation decay in the sheer 
dynamics; for the LLg it is very complicated [35]. Hence, while formally there might be a similar mechanism in gLe and sheer 
dynamics for generating subdiffusion, again, conceptually these dynamics are very different. 

To summarize this section, to our knowledge currently there is no stochastic model that fully reproduces the sheer diffusion. 
The superdiffusive sheer regime formally matches to diffusion known from Levy walks, as reproduced both by the LLg and 
CTRW theory, although in detail the parameter dependences for the MSD generated by both models are different. The subd¬ 
iffusive regime formally matches to what is generated by a gLe with power law correlated Gaussian noise. Conceptually all 
these stochastic models are very different from the sheer model, hence any similarity is purely formal and not supported from 
first principles. Based on this analysis it is tempting to conclude that one may need a correlated CTRW model to stochastically 
reproduce the sheer diffusion, or perhaps alternatively a simple Levy Markov chain model. 

The reason why we elaborated so explicitly on a possible stochastic modeling of the sheer dynamics is that our analysis 
might help to understand the reason for the controversy of how to stochastically model diffusion in polygonal billiards [3, 8, 
21, 22]. With the sheer map we are trying to capture a basic mechanism generating diffusion in these systems. However, as we 
showed above, the sheer diffusion seems to share features with generically completely different stochastic models depending 
on parameter variation. This might help to explain why different groups of researchers came to contradicting conclusions for 
modeling diffusion in polygonal billiards (not from first principles) by applying different types of stochastic processes. 


IV. CONCLUDING REMARKS 

In search of mathematically tractable deterministic models of normal and anomalous diffusion, which may shed light on the 
minimal mechanisms generating different transport regimes in non-chaotic systems, we have introduced a new model which we 
called the sheer map. This map mimics in a one dimensional space main features that distinguish periodic polygonal billiards 
from other models of transport, namely the complete absence of randomness and of positive Lyapunov exponents, and a sequence 
of splittings of a beam of particles due to collisions at singularities of the billiard walls. As observed in Refs. [12, 13, 23-26], 
in these cases the geometry determines the transport law, differently from standard hydrodynamics in which the geometry only 
yields the boundary conditions. Therefore the rule according to which the polygonal scatterers are distributed in space plays a 
crucial role. Here we have investigated the case of a specific deterministic rule modeling diffusion in polygonal billiards. 

In our one-dimensional sheer model S a the effect of the billiard geometry, which “slices” beams of particles ever more finely 
further and further away from the origin, is produced by the rate at which the size of the slices decreases with the position, 
i.e. by the value of a single control paraemeter a. For instance, a = 1/3 yields for the MSD and for the even higher order 
moments (A X 2 (n)) ~ n% and (AJ p (n)) ~ n^~ for long times. As we have discussed, the n 5 / 3 behavior coincides 
with the asymptotic MSD estimated numerically for a periodic polygonal channel made of parallel walls which form angles of 
90 degrees [12, 13] for which one would naively expect ballistic behaviour. 

It seems there do not exist too many models, neither in terms of deterministic nor stochastic dynamics, that exhibit sub-, super- 
and normal diffusive regimes under single parameter variation. The sheer map adds a new facet to this rather rare collection, 
as it generates all these different types of diffusion in a strictly deterministic and non-chaotic way. This suggests a mechanism 
explaining why in computer simulations of polygonal billiards so many different types of diffusion have been observed under 
parameter variation [12, 13, 23-26]. It may also help to explain the severe difficulties to model diffusion in such systems by a 
single, sufficiently simple anomalous stochastic process: We have argued that, depending on the value of its control parameter, 
the sheer diffusion matches mathematically to what is generated by very different classes of stochastic processes, which is in 
line with findings for polygonal billiard diffusion. 

It would be highly desirable to construct a simple stochastic process reproducing the full range of the sheer diffusion. It 
would also be important to extract a sheer map from a given polygonal billiard starting from first principles. This would enable 
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to check whether a slicing mechanism similar to the one proposed here in terms of a power law distribution of slicers is realistic. 
Correlation functions for the sheer map need to be calculated in order to fully appreciate its dynamics. And as the sheer map 
is analytically tractable, our model invites to play around with variations of the sheer idea for better understanding the origin of 
non-chaotic diffusion. 


Acknowledgements: 

The authors are grateful to Raffaella Burioni for illuminating remarks. LR acknowledges funding from the European Research 
Council, 7th Framework Programme (FP7), ERC Grant Agreement no. 202680. The EC is notresponsible for any use that might 
be made of the data appearing herein. CG acknowledges financial support from the MIUR through FIRB project “Stochastic 
processes in interacting particle systems: duality, metastability and their applications”, grant n. RBFR10N90W and the Institut 
Henry Poincare for hospitality during the trimester “Disordered systems, random spatial processes and some applications”. 


Bibliography 


[1] J.R. Dorfman. An introduction to chaos in nonequilibrium statistical mechanics. Cambridge University Press, Cambridge, 1999. 

[2] P. Gaspard. Chaos, scattering, and statistical mechanics. Cambridge University Press, Cambridge, 1998. 

[3] R. Klages. Microscopic chaos, fractals and transport in nonequilibrium statistical mechanics, volume 24 of Advanced Series in Nonlinear 
Dynamics. World Scientific, Singapore, 2007. 

[4] L. Rondoni and C. Mejia-Monasterio. Nonlinearity, 20:R1, 2007. 

[5] P. Castiglione, M. Falcioni, A. Lesne, and A. Vulpiani. Chaos and Coarse Graining in Statistical Mechanics. Cambridge University 
Press, Cambridge, 2008. 

[6] U. Marini Bettolo Marconi, A. Puglisi, L. Rondoni, A. Vulpiani. Phys. Rep., 461:11, 2008. 

[7] O.G. Jepps, L. Rondoni. J. of Phys. A, 43:P05015, 2010 

[8] G.M. Zaslavsky. Phys. Rep., 371:461, 2002. 

[9] R. Klages, G. Radons, and I.M. Sokolov, editors. Anomalous Transport: Foundations and Applications. Wiley-VCH, Berlin, 2008. 

[10] O.G. Jepps, S.K. Bathia, D.J. Searles. Phys. Rev. Lett. 91:126102, 2003 

[11] A. Igarashi, L. Rondoni, A. Botrugno, M. Pizzi. Commun. Theor. Phys. 56:352, 2011 

[12] O.G. Jepps and L. Rondoni. J. Phys. A: Math. Gen., 39:1311, 2006. 

[13] O.G. Jepps, C. Bianca, and L. Rondoni. Chaos, 13:013127, 2008. 

[14] I.M. Sokolov. Soft Matter 8:9043, 2012. 

[15] R. Metzler and J. Klafter. Phys. Rep., 339:1, 2000. 

[16] G. Zumofen and J. Klafter. Phys. Rev. E, 51:1818, 1995. 

[17] N. Korabel and R. Klages. Phys. Rev. Lett., 89:214102, 2002. 

[18] N. Korabel and R. Klages. Physica D, 187:66, 2004. 

[19] T. Manos and M. Robnik. Phys. Rev. E, 89:022905, 2014. 

[20] B. Li and L. Wang and B. Hu. Phys. Rev. Lett., 88:223901, 2002 

[21] B. Li, J. Wang, L. Wang, and G. Zhang. Chaos, 15:015121, 2005. 

[22] Carl P. Dettmann. Commun. Theor. Phys., 62:521, 2014. 

[23] R. Artuso, I. Guarneri, and L. Rebuzzini. Chaos, 10:189, 2000. 

[24] D. Alonso, A. Ruiz, and I. de Vega. Physica D, 187:184, 2004. 

[25] M. Schmiedeberg and H. Stark. Phys. Rev. E, 73:031113/1-9, 2006. 

[26] D.P. Sanders and H. Larralde. Phys. Rev. E, 73:026205, 2006. 

[27] J. Hardy and J. Weber. J. Math. Phys., 21:1802, 1980. 

[28] V. Delecroix, P. Hubert, and S. Lelievre. Ann. Scient. Ec. Norm. Sup. 4e serie, t. 47, 1085 (2014). 

[29] J.H. Hannay and R.J. McCraw. J. Phys. A: Math. Gen., 23:887, 1990. 

[30] E. Gutkin. J. Stat. Phys., 83:7, 1996. 

[31] G.M. Zaslavsky and M. Edelman. Chaos, 11:295, 2001. 

[32] C.P. Dettmann and E.G.D. Cohen. J. Stat. Phys., 101:775, 2000. 

[33] F. Cecconi, M. Falcioni, A. Vulpiani, and D. del Castillo-Negrete. Physica D, 180:129, 2003. 

[34] S. Denisov, J. Klafter, and M. Urbakh. Phys. Rev. Lett., 91:194301, 2003. 

[35] E. Barkai and V. Fleurov. J. Stat. Phys., 96:325, 1999. 

[36] E. Barkai, V. Fleurov, and J. Klafter. Phys. Rev. E, 61:1164, 2000. 

[37] R. Burioni, L. Caniparoli, and A. Vezzani. Phys. Rev. E, 81:060101(R), 2010. 

[38] P. Barthelemy, J. Bertolotti, and D.S. Wiersma. Nature, 453:495, 2008. 

[39] N. Kumar, U. Harbola, and K. Lindenberg. Phys. Rev. E, 82:021101, 2010. 

[40] U. Harbola, N. Kumar, and K. Lindenberg. Phys. Rev. E, 90:022136, 2014. 


14 



[41] J.M. Porra, K.-G. Wang, and J.Masoliver. Phys. Rev. E, 53:5872, 1996. 

[42] A. Taloni, A.V. Chechkin, and J. Klafter. Phys. Rev. Lett., 104:160602, 2010. 

[43] J.-H. Jeon and R. Metzler. Phys. Rev. E, 81:021103, 2010. 

[44] A. Rebenshtok, S. Denisov, P. Hanggi, and E. Barkai. Phys. Rev. Lett., 112:110601, 2014. 

[45] A.V. Chechkin, F. Lenz, and R. Klages. J. Stat. Mech.: Theor. Exp., 2012:L11001, 2012. 

[46] The slicer odd moments vanish because of the symmetry of the map considered here. If the moments are computed only over the positive 
half of the chain, also the odd moments can be identified with those of [37], 


15 


